Predicting suitable habitat for the endangered tree Ormosia microphylla in China

Climate change has significantly influenced the growth and distribution of plant species, particularly those with a narrow ecological niche. Understanding climate change impacts on the distribution and spatial pattern of endangered species can improve conservation strategies. The MaxEnt model is widely applied to predict species distribution and environmental tolerance based on occurrence data. This study investigated the suitable habitats of the endangered Ormosia microphylla in China and evaluated the importance of bioclimatic factors in shaping its distribution. Occurrence data and environmental variables were gleaned to construct the MaxEnt model, and the resulting suitable habitat maps were evaluated for accuracy. The results showed that the MaxEnt model had an excellent simulation quality (AUC = 0.962). The major environmental factors predicting the current distribution of O. microphylla were the mean diurnal range (bio2) and precipitation of the driest month (bio14). The current core potential distribution areas were concentrated in Guangxi, Fujian, Guizhou, Guangdong, and Hunan provinces in south China, demonstrating significant differences in their distribution areas. Our findings contribute to developing effective conservation and management measures for O. microphylla, addressing the critical need for reliable prediction of unfavorable impacts on the potential suitable habitats of the endangered species.


Establishing species occurrence records
Occurrence data of O. microphylla in China were obtained from the Chinese Virtual Herbarium (http:// www.cvh.ac.cn), Plant Photo Bank of China (http:// www.Plant photo photo.cn), National Specimen Information Infrastructure (http:// www.nsii.org.cn/), field surveys, and the published literature.Each occurrence record was scrutinized for data quality and suitability with reference to the model's requirements.Records with unclear latitude and longitude information and duplicated distribution points were discarded.A total of 45 verified distribution points of O. microphylla were kept for final analysis to generate the potential distribution based on MaxEnt modeling.The occurrence-point data were stored in a CSV file, sorted by species name, longitude and latitude, and plotted on a map (Fig. 1).

Selecting environmental variables
Data on 19 bioclimatic parameters, including monthly temperature and precipitation, were selected and downloaded from the WorldClim website (http:// www.world clim.org) (Table 1).The comprehensive climatic data from 1950 to 2000 were based on monthly meteorological data harvested from different weather stations worldwide.The raw climatic data were interpolated according to a grid with a 2.5-arcminute resolution to generate the global spatial dataset.The rather general annual temperature and precipitation variables were deleted from the dataset 23,24 .Strong correlation among the environmental variables can lower model accuracy.Therefore, Pearson's correlation coefficients (r) were calculated between variable pairs.Variables reaching |>0.8| were removed to avoid multicollinearity impacts 25 .Ultimately, we shortlisted eight variables to construct the MaxEnt model for O. microphylla in this study.
To avoid model overfitting, the "ENMeval" software package in R 4.0.2 was employed to optimize MaxEnt model 26 .For optimization, regular Multiplier (RM) and Feature Class (FC) were included in MaxEnt.RM could smooth the model and minimize model over-fitting.FC, corresponding to the response type of suitability values to each variable, could determine the potential shape of response curves, including Linear (L), Product (P), Quadratic (Q), Threshold (T), and Hinge (H) 27 .In the process of parameter optimization, the RM value ranged from 0.5 to 4.0 with an increment of 0.5.Six combinations of FCs (L, LQ, LQP, LQT, QPT, PHT) were selected.Combined with the 8 values of RM, this resulted in a total of 48 parameter combinations.The 48 parameter combinations were input into ENMeval for comprehensive testing.We used the delta from the Akaike information criterion (AIC), AICc, the difference between training and testing AUC (AUC.DIFF) and the 10% training omission rate (OR10) to evaluate the model's fitting degree and complexity on species distribution 28 .

Modeling species distribution
The eight selected environmental variables and species occurrence records of O. microphylla were loaded into MaxEnt 3.3.Then, 75% of the distribution points were randomly selected as testing data to establish the prediction model, and the remaining 25% were used to verify the model's accuracy.To avoid the instability caused by randomly selected data, the initial operation of the model was repeated ten times.A Jackknife test was applied to measure the relative importance of each variable on the distribution of this species by percent contribution to the overall model fit and permutation importance.Permutation importance represents how heavily the final model depended on a certain variable 29 and was calculated by random values of each variable between presence and background points and then measures resulting drop in AUC with a larger drop indicating higher importance of that variable to the model 30,31 .Then, the receiver operating characteristic (ROC) curve and the area enclosed by the abscissa as the area under the curve (AUC) were used to evaluate model accuracy 24 .The AUC value ranges  from 0.5 for an uninformative model to 1 for perfect discrimination 32 .An AUC value close to 1 indicates more habitat suitability deviation from a random value and a stronger correlation between environmental variables and habitat classification, signifying a better model prediction 33,34 .The suitable habitat maps were calculated using the logistic output of MaxEnt, which ranges from 0 to 1.The outputs of the MaxEnt model in logistic format were imported into QGIS software for data extraction.China's map served as the base to extract the potential suitable distribution of O. microphylla.The average logical value combined with the actual distribution were used to classify distribution-level values and corresponding distribution ranges.The specific suitability is divided into five categories: fail (0-0.15),poor (0.15-0.3), fair (0.3-0.45), good (0.45-0.6), and excellent (>0.6) 35,36 .Different suitable levels indicate the differential potential distribution probability of the species in specific areas.The higher the probability, the more suitable is an area for species occurrence.

Evaluating model performance
The current distribution areas of O. microphylla were predicted using the MaxEnt model (Fig. 2).The AUC value for the MaxEnt models was 0.962, which was significantly higher than the AUC value of a random prediction (0.5).This result indicated that the prediction results were "excellent", and the MaxEnt model was reliable for predicting the potential geographical distribution areas of O. microphylla in China.

Key environmental factors and validation of modeling results
The key environmental factors shaping the potential distribution of O. microphylla were determined according to their contributions to the MaxEnt modeling process by the jackknife test (Table 2).It found that bio14 (precipitation of the driest month), bio4 (temperature seasonality), and bio16 (precipitation of the wettest quarter) had the highest contribution of 53.3%, 22.6% and 12.2%, respectively (aggregate contribution of 88.1% to the  The results of jackknife tests indicated that when only individual variables were used, the regularized training gain, test gain, and AUC values of the mean diurnal range (bio2) and precipitation of the driest month (bio14) declined the most.These two variables, demonstrating the strongest influence on species distribution, were the main factors influencing the geographical distribution of O. microphylla (Fig. 3).
The species response curve depicts the relationship between environmental variables and the probability of species incidence.They show the target species' biological tolerances and habitat preferences.The distribution probability of O. microphylla increased with the value of each environmental variable within a certain range.It decreased with the increase of the variable after reaching a certain peak value (Fig. 4).Based on the species response curves, O. microphylla prefers the mean diurnal range (bio2) range from 7 to 8°C, and precipitation of the driest month (bio14) ranges from 30 to 40 mm (Fig. 4).

Predicting suitable habitats of O. microphylla in China
The five categories of O. microphylla suitability habitats for the current climate were mapped (Table 3 and Fig. 5).The suitable areas for the fail, poor, fair, good, and excellent suitability categories occupied 183.67×10 4 km 2 , 41.76×10 4 km 2 , 23.19×10 4 km 2 , 17.42×10 4 km 2 , and 16.55×10 4 km 2 , respectively, comprising 19.13%, 4.35%, 2.42%, 1.81%, and 1.72% of China's total land area (about 960×10 4 km 2 ) (Table 3).The excellent suitable habitat areas were mainly distributed in the south China provinces of Guangxi, Guizhou, Hunan, Guangdong, and Fujian, covering approximately 14.59×10 4 km 2 (Table 3).However, the provinces contiguous or proximal to the main distribution areas, namely Hubei, Zhejiang, and Taiwan, have no distribution records in the good and excellent categories.This pronounced discrepancy indicates that the potential range of O. microphylla is considerably larger than the actual distribution range, and the species has the potential to expand the current range to fill the potential areas.

Discussion
This is the first study employing the MaxEnt model in conjunction with QGIS to analyze the spatial distribution pattern of O. microphylla and climate change impacts on its potential geographical range, and predict its potential distribution area in China quantitatively and intuitively.We built the model with 45 valid occurrence records and eight key environmental factors (Fig. 1).The AUC value generated by MaxEnt reached 0.962 (Fig. 2), indicating the prediction had a high degree of fit between the climatic variables and the actual suitable habitat of O. microphylla.The model jackknife test provided a measure of the contribution of the various environmental variables in influencing the suitability for O. microphylla.The prediction results provide a theoretical basis for developing protection and management measures for the endangered species.This study analyzed the relationship between the occurrence probability of O. microphylla and key environmental variables.The results showed that precipitation of the driest month (bio14) was the primary bioclimatic variable affecting the presence of O. microphylla, with optimum conditions at 30 to 40 mm.As a plant growth prerequisite, precipitation is the primary limiting factor for almost all species 37,38 .Variations in precipitation and consequent changes in temperature and humidity disrupt the balance of soil moisture and most physiological plant functions 36,39 .The driest month in China occurs in March, April, or May, accompanied by higher temperatures.Moreover, this Spring season is critical for the germination of O. microphylla seeds.Sufficient precipitation  2 and Fig. 3), with optimum conditions at 7-8 °C (Fig. 4).The mean diurnal range denotes the temperature variations in a day, which may impact plant growth, especially daytime photosynthesis and respiration that contribute to nutrient buildup 40 .Some pertinent processes associated with photosynthesis are highly sensitive to temperature.Low temperatures can reduce the hydrolysis and transport of starch accumulated within the chloroplast to depress the photosynthetic rate.It can also suppress the germination rate and biomass accumulation and modify the survival strategies.On the contrary, the photosynthetic rate increases in response  to temperature rise until reaching a thermal optimum, after which the rate declines due to enzyme deactivation at high temperature 41,42 .Other studies have shown faster temperature rises at nighttime than daytime, which is not conducive to nutrient accumulation 43,44 .However, the specific mechanism concerning the influence of mean diurnal range on O. microphylla is unclear.Therefore, research on this topic could be strengthened in the future.The other variables, such as bio4 (temperature seasonality) and bio16 (precipitation of the wettest quarter), have an important influence on the habitat suitability distribution of O. microphylla (Table 2).Temperature seasonality represents the temperature variations in a year.A larger standard deviation brings a greater coefficient of variation 45 .Therefore, this species is not suitable for an environment with a wide annual temperature amplitude.The wettest period in China occurs in June, July, and August, which also have the highest yearly temperatures 46 .Moreover, this Summer season is also a critical period for the maturation of O. microphylla fruits and seeds.Sufficient precipitation in the wettest quarter can raise atmospheric humidity and soil moisture to foster fruit growth and maturation.The moist conditions also increase fats and soluble sugars in seeds, boosting the energy supply for germination and initial seedling growth.The ability of seedlings in different habitats to resist harsh environmental conditions will be correspondingly enhanced 47 .
Under current climate conditions, the excellent suitable habitats for O. microphylla mainly cover Guangxi, Guizhou, Hunan, Guangdong, and Fujian provinces in south China (Table 3 and Fig. 5).The model's prediction matches the existing data, highlighting the accuracy of the prediction.However, the predicted suitable habitat is much larger than the current known species distribution.No occurrence record is found in Sichuan, Yunnan, Hubei, Anhui, Zhejiang, Jiangsu, and Taiwan provinces, even though they are the species' potentially important habitats.This notable discrepancy may be related to limited research on this species, resulting in inadequate data to describe species environmental requirements accurately.On the other hand, MaxEnt has the inherent trait of evaluating only niche-based species presence data 48 .The model predicts the species' fundamental niche rather than the actual niche, resulting in the predicted potential distribution area being larger than the actual range.However, the species may fail to disperse to some suitable areas due to biogeographical barriers, such as human interference, topographical obstacles, and inter-specific competition 49 .
Based on the predicted results, the suitable habitat of O. microphylla is relatively widespread.However, in reality, the prospect of population development is not optimistic.To prevent the extinction of this endangered plant, in situ conservation of the known O. microphylla community is necessary to minimize the adverse effects of human activities on species survival.More distribution point information can be obtained by conducting O. microphylla studies to deepen the knowledge of biological characteristics, genetic structure, and artificial cultivation methods, as well as investigating and monitoring the apparent "omission" areas.Collecting seeds for artificial breeding and selecting potential suitable areas with less human activity for wild release can expand the spread and size of wild populations.
The ability of the MaxEnt model to use scarce distribution records to infer the environmental tolerance and ecological niches of O. microphylla is supported by the generation of highly accurate predictions.However, it was difficult for the species to establish and survive in all suitable areas.The distribution areas of its communities have been progressively reduced despite the expansion of suitable habitats 50 .Besides external environmental and climatic conditions, other factors may have important implications for species distribution at different spatial scales 51 .They include physiological constraints, topography, soil, human disturbance, spatial constraint, dispersal mode, competition, and response to external factors and drivers.For instance, the seeds of O. microphylla are hard and compact, with a less permeable cortex and a long dormancy period.Thus, its germination rate is extremely low, making it difficult to survive in suitable areas 21 .In addition, land use and land cover changes caused by human activities have significantly impacted species distribution 52 .However, these factors take time to take effect, which may cause species distribution to lag behind climate change 53 .
Under the climatic conditions of increasing CO 2 concentration and temperature in the future, the annual mean precipitation may increase across China 54 .Our research showed that the main factors influencing the distribution of O. microphylla were the mean diurnal range and precipitation of the driest month.With the increase of temperature and precipitation in the future climate change, sufficient soil moisture and suitable soil temperature can better meet the needs of seed germination and vegetation growth 55 .Therefore, the suitable area of O. microphylla may vary with the magnitude of future climate change.
In short, predicted results based on the MaxEnt model could help identify additional localities where the O. microphylla may already exist but have not yet been detected.Alternatively, the predicted suitable areas offer new habitats to expand the species range or are used as priority areas for introducing and cultivating this rare tree.Therefore, deeper cognate research should be conducted to expand the knowledge and practice base.The locations of key protected areas for O. microphylla could be identified to strengthen their supervision and management and minimize human activities.Uninformed and subjective designation of protection areas can be avoided.The numerical range of the main environmental variables can provide an objective reference to select artificial breeding and rescue sites.Thus, the protected plants can be nurtured in appropriate environmental conditions to rationalize resource development and utilization.Moreover, this study indicates that other aspects can be improved.In addition to the selected bioclimatic variables, other factors such as soil, water, land use change, different climate scenarios, and dispersal ability may affect the prediction results.Therefore, future studies can comprehensively consider a broader range of environmental factors to develop a response strategy for O. microphylla to tackle climate change impacts.

Conclusion
The suitable areas of O. microphylla in China have been accurately predicted by MaxEnt modeling using mainly bioclimatic variables.The excellent and good suitable habitats were found primarily in China's southern provinces, including Guizhou, Hunan, Guangxi, Guangdong, and Fujian provinces.The crucial environmental

Figure 1 .
Figure 1.The locations of 45 verified current occurrence points of O. microphylla (white circles) in central and south China with reference to elevation.The map was prepared by Lijuan Wei, Guohai Wang and Chunping Xie in QGIS 3.34.0(https:// www.qgis.org/ en/ site/).

Figure 2 .
Figure 2. The validation of the MaxEnt model predicting O. microphylla distribution: (a) Omission rate; and (b) ROC curve.

Figure 3 .
Figure 3.The jackknife test for evaluating the relative importance of environmental variables for O. microphylla.

Figure 4 .
Figure 4. Response curves of the environmental variables to distribution probability.

Figure 5 .
Figure 5.The predicted O. microphylla distribution range classified into five suitability categories obtained by MaxEnt modeling.The map was prepared by Lijuan Wei, Guohai Wang and Chunping Xie in QGIS 3.34.0(https:// www.qgis.org/ en/ site/).

Table 1 .
List of 19 environmental variables considered at the initial stage of model development.Eight variables with the code shown in bold font were chosen for the MaxEnt modeling study.

Table 2 .
Percent contribution and permutation importance levels of the eight environmental variables included in the MaxEnt models, ranked by percentage contribution.

Table 2 )
. The mean diurnal range (bio2, 45.8%) and precipitation of the driest month (bio14, 38%) had the highest permutation importance.The correlation coefficients (r) between the eight environmental factors were below 0.8.The two bioclimatic factors (bio14, bio2) were identified as the main drivers of the modern geographical distribution of O. microphylla.

Table 3 .
Predicted suitable areas for O. microphylla under the current climate scenario in various provinces (10 4 km 2 ).The ratio denotes the predicted suitable area divided by the total land area of the respective province or autonomous.